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We investigate the stability of neutrino mass matrix predictions on important and currently 
unknown observables. Those are the octant of 623 , the sign of sin <5 and the neutrino mass ordering. 
Determining those unknowns is expected to be useful in order to distinguish neutrino mass models. 
Therefore it may be interesting to know how robust the predictions of a mass matrix for the octant of 
023 or the neutrino mass ordering are. By applying general multiplicative perturbations we explicitly 
quantify how probable it is that a perturbed mass matrix predicts an octant of 023 different from 
the original mass matrix, or even a neutrino mass ordering different from the original one. Both 
the general case and an explicit flavor symmetry model are studied. We give the probabilities as 
a function of the smallest neutrino mass, showing that for values exceeding 0.1 eV the chance to 
switch the prediction quickly approaches 50 %. 
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I. INTRODUCTION 

In recent years a consistent picture of lepton mixing has 
emerged with several parameters being determined 
with increasing precision (for a recent global fit of all 
existing data, see Ref. [2]). A remarkable pattern has 
emerged, with one close-to-maximal mixing angle, one 
large and one small mixing angle, the latter being of the 
order of the largest quark mixing angle. While the over¬ 
all picture of the leptonic mixing matrix is clear, com¬ 
parable precision with respect to the quark sector is still 
lacking, but future experiments and facilities exist that 
will improve the errors on the parameters by remarkable 
amounts, see e.g. [3]. Of particular interest in neutrino 
physics are the octant of the atmospheric neutrino mix¬ 
ing angle 023 and of course sin S, the parameter governing 
leptonic CP violation. The mass ordering and the value 
of the smallest neutrino mass are also unknown (while 
not yet determined, we will assume here that neutrinos 
are Majorana particles). 

The astonishing disparity between lepton and quark 
mixing has lead to huge efforts in flavor symmetry model 
building |M|. Many neutrino mixing schemes have been 
proposed (see e.g. H), and many models exist that can 
generate these schemes. The question is now of course to 
distinguish the various models or scenarios and identify 
the correct one. One could expect that the determina¬ 
tion of the unknown neutrino parameters, in particular 
the sign of sin 5, the octant of 023 or the mass ordering 
will be crucial. In this paper we analyze how robust these 
parameters are with respect to perturbations of the mass 
matrix. Perturbations to a mass matrix are expected to 
be present because of various reasons, e.g. renormaliza¬ 
tion effects including thresholds, misalignment of the vac¬ 
uum expectation values of the ffavons which are crucial 
in flavor symmetry models, non-canonical kinetic terms, 
higher-dimensional operators, etc. By quantifying how 
probable it is that a perturbed mass matrix changes its 
predictions for a currently unknown neutrino parameter, 
one can estimate how robust the predictions are. Analyz¬ 


ing this issue is the purpose of the present paper. As the 
probability to change the predictions depends strongly on 
the neutrino mass scale, this is especially crucial for siz¬ 
able neutrino masses, with the extreme case being quasi¬ 
degenerate neutrino masses. 


Our procedure is as follows: we start with a large set 
of mass matrices that are allowed according to current 
global fits, but have a certain property that is of in¬ 
terest to us, say, 023 < 7r/4. Then we multiplicatively 
perturb the mass matrices in a general way, and check 
how many percent of the resulting mass matrices change 
the property of interest, i.e. predict 023 > 7r/4 after per¬ 
turbation. This percentage is a function of the smallest 
neutrino mass. We demonstrate that, as one may ex¬ 
pect, the percentage increases strongly with the smallest 
neutrino mass and is in general larger for the inverted 
ordering than for the normal one. The solar neutrino 
mixing angle is subject to the largest instability among 
the mixing angles, the CP phase 6 as well. The sign of 
sin (5 is more likely to change than the octant of 623 - For 
values of the smallest neutrino mass around 0.1 eV and 
larger, even the mass ordering can change from normal 
to inverted. In general, for neutrino masses larger than 
0.1 eV the probability to change a prediction quickly ap¬ 
proaches 50%. While intuitively many findings are ex¬ 
pected, there has never been a quantitative study ad¬ 
dressing these issues. We also analyze an explicit model 
based on A 4 , in which correlations among the mass ma¬ 
trix parameters are present. We And qualitatively similar 
results. This demonstrates the challenge to distinguish 
neutrino mass models or scenarios unless corrections are 
taken into account. 


The paper is build up as follows: in Section |TT] we 
present the procedure and then discuss the perturbation 
of a general mass matrix. Section [ni| deals with a specific 
A 4 model, before we conclude in Section EYi 
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II. PERTURBATION OF A GENERAL MASS 
MATRIX 

A. Method 

Let us start with a zeroth-order neutrino mass matrix 
Mq, constructed by 

Mo = C/odiag(mi,TO 2 ,m 3 ) , 

where Ug, parametrized as usual [T], includes the Majo- 
rana phases and (mi, m 2 , m 3 ) can be determined by the 
mass-squared differences (following the definition of Ref. 
[2]) 5m^ = 1712 — ml, Am^ = m| — (mf -I- ml)f2. The 
lightest neutrino mass is mi for the normal mass order¬ 
ing, m 3 for the inverted one. 

We consider now a general multiplicative perturbation 
to the individual mass matrix entries: 

(Mo)a^ —>■ Map = {Mo)ap (1 + £ap) , (1) 

where tap = £pa are six small complex numbers. Note 
that multiplicative perturbations are conservative, one 
could also add terms eap M™“ to each entry, i.e. cor¬ 
rections proportional to the largest entry in the mass 
matrix. Such additive corrections are expected to give 
qualitatively similar perturbations as the ones we will de¬ 
rive here. However, they will be at least as sizable as the 
multiplicative ones under study, as they influence small 
entries of the mass matrix more significantly (note that 
with multiplicative corrections texture zeros and the as¬ 
sociated correlations they introduce are not significantly 
changed). In addition, often and extensively studied cor¬ 
rections from the charged lepton sector could be included 
as well. We have nothing new to add to this aspect, and 
in addition those correction are model-dependent, and 
furthermore independent on neutrino mass and ordering. 
We stick in the present paper to the conservative case of 
multiplicative corrections to mass matrices and the anal¬ 
ysis thereof. 

One has several possibilities to choose the initial pa¬ 
rameters |H]. We decided to choose in Mq the mixing 
angles and mass-squared differences ( 0 i 3 , 0 J 21 ^23 

(Am^)°) randomly within their current 3cr con¬ 
fidence intervals, while for both Dirac and Majorana 
CP phases, we randomly generate them in [0,27r]. We 
will however be interested in a certain property, say 
6*23 < 7 r/ 4 . Therefore, this condition is imposed on Mq. 
After Mq is constructed, we randomly generate the six 
tap with ^ |e| < 0.2 and ^ |ep > 0.01^. We require that 
M after perturbation (having mixing angles @13, 612, 823 
and mass-squared differences Sm^, /S.m?) is still compat¬ 
ible with current data within 3a. We are interested in 
the percentage of successful neutrino mass matrices M 
that are within 3a, but have went from 023 < 7 r /4 to 
023 > 7 r/ 4 . Put another way, we obtain the probability 
for the perturbed mass matrix to change the characteris¬ 
tic prediction we are interested in. The same procedure 
is performed for the sign of sinJ and for the mass or¬ 
dering. We are interested in how the results depend on 


the smallest neutrino mass. To make robust statements, 
we want 10000 successful mass matrices for each value of 
the smallest mass. Hence, the numerical analysis is quite 
CPU-intensive in particular for neutrino masses near and 
above 0.1 eV. 

A comment on the upper and lower limits on the e is in 
order: a compromise between a “reasonable” percentage 
of successful mass matrices after perturbation on the one 
hand, and guaranteed perturbations to the mixing pa¬ 
rameters on the other hand, needs to be found. A lower 
limit on the perturbations is needed because if we have no 
lower limit the vast majority of successful mass matrices 
is essentially identical to the original ones. Larger upper 
limits on the eap than the ones we use increase the CPU¬ 
time for sizable neutrino mass significantly. The limits 
on the Cap might be interpreted similarly to the model 
analysis in Section |HI[ namely as VEV misalignment in 
a flavor symmetry model of order a few percent. With 
typically 3 to 5 VEVs playing a role, see Eqs. ([^[^, the 
upper limit on the sum of X) kl ^ could be under¬ 
stood. We would like to avoid too much cancellations in 
the e, hence a lower limit should be present. Threshold 
effects with RG corrections might be another interpreta¬ 
tion of the e. We prefer however to stay here as model- 
independent as possible. In any interpretation of the e, 
a given model might induce a correlation between them. 
This is realized in the model that is studied in Sec. mil 

Anyway, we have checked that for small neutrino 
masses, where the analysis takes still reasonable CPU¬ 
time, the results do hardly depend on the precise values 
of the upper and lower limits of the e, up to longer CPU¬ 
time for larger upper limits. With this check we gained 
confidence in our choice of limits. 

A few words on generating the events: the most 
straightforward method to realize it is simply "generate- 
and-reject", which means to generate enough events 
without the constraint and then reject those which vi¬ 
olate our constraints. This is of low efficiency especially 
for large values of the smallest neutrino mass, since M af¬ 
ter perturbation is quite likely to go out of the 3cr bound 
for a quasi-degenerate mass spectrum. Therefore, besides 
optimization of the algorithm which includes fast diago- 
nalization of M and extracting the neutrino parameters, 
we use the "generate-and-tune" method: we first ran¬ 
domly produce six e and then choose their phases such 
that the following y^-function is minimized: 



where pi are the oscillation parameters which are irrel¬ 
evant. For example when studying the stability of 823 , 
we take 5m^, Am^, 0 i 3 and 0 i 2 as irrelevant parameters; 

and ai denote respectively the best-fit values and la 
errors of the corresponding parameters from Ref. |5]. We 
have checked that the events generated in this way have 
almost the same distribution as those generated by the 
"generate-and-reject" method, but the procedure is more 
efficient and faster. 
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One comment should be added here: often the mass 
matrices that are resulting from flavor symmetry models 
have a feature called "form-invariance", i.e. the eigen¬ 
values do not depend on the mixing angles (infamous 
tri-bimaximal mixing is one particular example) and our 
analysis might be irrelevant in this case. However, if 
breaking terms are added to the mass matrices the form- 
invariance is lost. 


B. Results 


We first look at the correlation of the mixing angles 
by simply plotting 023)^13)^12 and 6 after perturbation 
against the original mixing angles 023) ^i3> ^12 and 5 ° be¬ 
fore perturbation. This should give us as feeling on how 
much the perturbation changes the mixing angles. One 
expects that 0i2, being related to the smaller of the two 
mass-squared differences, will be most unstable. One also 
suspects S, being related to phases of various mass ma¬ 
trix elements, to be quite unstable. Furthermore, the 
larger the smallest neutrino mass, the larger the average 
perturbation. 

The result is shown in Figs. min] and (to illus¬ 
trate the outcome in a optimal way, we use an upper 
bound kl < 0-04 instead of 0 . 2 ). The plots confirm 
the expectation. When the smallest mass is 0.001 eV, 
023 is stable, typically deviating from its original value 
by ~ 1 ° (depending on the upper bound of ^ |e|). When 
the smallest mass increases, the points spread and for 0.1 
eV there is hardly any correlation left. This conclusion 
equally applies for 0i3, as shown in Fig. Note that it is 
the most precisely measured angle, and the range of the 
y-axis is much narrower than for 023 . However, 0i2 and <5 
are very unstable even for small masses, as can be seen in 
Figs. 1^ and 0 (note the large range of the axes for the plot 
with 5 ). This implies that distinguishing models based 
on precision measurements of 0i2 and/or S is not a par¬ 
ticularly reliable method unless corrections are carefully 
included in the predictions of a model. Recall that the 
plots are for the normal mass ordering. For the inverted 
ordering, d and 0i2 will be uncorrelated with 5 *^ and 032 
even for the smallest value of 0.001 eV, while 023.13 are 
slightly more uncorrelated (see below). 

After these preliminaries, we evaluate now the percent¬ 
ages of the perturbed mass matrices that change the oc¬ 
tant of 023, the sign of sini 5 or the mass ordering (we 
do not consider Majorana phases, as their experimental 
determination is questionable). The results are shown in 
Figs. EE and In the plots we indicate two interest¬ 
ing mass scales \Jbrn? ~ 0.008 eV and VAm^ ~ 0.05 eV, 
which will be discussed in detail later. We also plot the 
rather strong 95 % C.L. limit on neutrino masses (com¬ 
bining various cosmological data sets) as given by Planck 
< 0 . 23 eV. 
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Figure 1. Correlation of ^23 (after perturbation) with 623 
(before perturbation). The lightest neutrino mass is (top to 
bottom) 0.001, 0.04 and 0.1 eV. 



Figure 2. Correlation of 0 i 3 (after perturbation) with 0 J 3 
(before perturbation). The lightest neutrino mass is (top to 
bottom) 0.001, 0.04 and 0.1 eV. 
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Figure 3. Correlation of 0 i2 (after perturbation) with 0i2 
(before perturbation). The lightest neutrino mass is (top to 
bottom) 0.001, 0.04 and 0.1 eV. 



Figure 4. Correlation of sin 5 (after perturbation) with sini?'^ 
(before perturbation). The lightest neutrino mass is (top to 
bottom) 0.001, 0.04 and 0.1 eV. 
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Figure 5. Percentage of perturbed mass matrices that give 
023 in the second octant when the unperturbed mass matri¬ 
ces give 023 in the first octant (blue for normal mass order¬ 
ing, red for inverted). The light-gray and gray regions start 
at V 5m? ~ 0.008 eV and \/Am^ ~ 0.05 eV, respectively, 
while the yellow region represents the strongest cosmology 
constraint on neutrino masses [9]. 



smallest mass/eV 


Figure 6. Percentage of perturbed mass matrices that give 
negative sin 5 when the unperturbed mass matrices give posi¬ 
tive sin 5 (blue for normal mass ordering and red for inverted). 
For other details, see Fig. 



Figure 7. Percentage of perturbed mass matrices that give 
the inverted or normal mass ordering when the unperturbed 
mass matrices correspond to the normal (blue) or inverted 
(red) mass ordering, respectively. For other details, see Fig. 
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C. Discussion 

Very simple arguments are enough to understand the 
features of the results. 

In Figs. mil!] we have indicated two relevant neu¬ 
trino mass scales, \/brn? ~ 0.008 eV and VAm^ ~ 
0.05 eV. Let us consid er th e normal mass ordering. Be¬ 
low the mass scale 'JSw? ~ 0.008 eV, all three mix¬ 
ing angles and the CP phase should be rather stable: 
neither 6 m?/{m\ + m\), associated to 12 -mixing, nor 
-I- TO§), associated to 13- and 23-mixing, are 
small. However, when the smallest mass increases, first 
6 m?/{m\ -I- m\) and then -I- m\) decrease and 

become small. Consequently first 6*12 and then 6 * 13^23 
will become unstable. Increasing the smallest mass fur¬ 
ther, corresponding to more and more quasi-degenerate 
masses, leads to all mixing angles becoming very unstable 
under perturbations. 

From Fig. we can see that for a smallest mass of 
0.04 eV, 012 is unstable (since 6 m?/{m\ + m|) is small) 
as there is no significant correlation between 0 i 2 and 
012 . In contrast to 0 i 2 , 0 i 3 and 023 are more stable, as 
shown in Figs. and Increasing the smallest mass, 
023 and 013 become unstable when Ato^/(toi -I- m|) be¬ 
comes sm all, w hich happens when the smallest mass goes 
beyond VAm^ ~ 0.05 eV. 

More quantitative is Fig. The probability of 

changing the octant is about 10 % for s mall m asses and 
remains constant until the mass scale y/Am? ~ 0.05 eV 
is reached. Approaching quasi-degenerate masses gives 
a probability of almost 50 % to change the octant, i.e. 
the octant is random and maximally unstable. 

The predictions are more stable for the normal mass 
ordering than for the inverted one. This conclusion 
comes from the comparison of the relevant percentages 
in Figs. and Note that for a vanishing smallest 
mass we have (toi, mf,m.|) « {Q, 6 m?,Am?) for a 

normal ordering whereas for the inverted case we have 
(TOi,m|,TO§) « {Am?, Am? + 6m‘^,0). Therefore, 
6 m? I{m?[ -fm^) is small from the beginning and of order 
6 rr? jAm? ~ 0.03. Indeed the probability to change 
the octant starts with about 15% and increases when 
Am? /{m?i -I- to|) becomes small for smallest masses 
of 0.05 eV and larger. Obviously for quasi-degenerate 
neutrino masses there will be no difference between the 
mass orderings. 

Comparing Fig. and we can see that between 
0.01 eV to 0.1 eV, the probability of sin <5 changing 
its sign is larger than the probability of 023 changing 
its octant. This is caused by the fact that phases of 
eigenvectors in a diagonalization procedure are always 
more sensitive to perturbation than their absolute 
values. One might also argue that 6 is related to the 
Jarlskog invariant J = Im(t7ei U *2 U{j,i) which is 
proportional to sin 5 sin 0 i 2 sin 0 i 3 sin 023 , which means 
that (for normal ordering) the probability of sin 6 


changing its sign should be similar to the probability 
of 023 changing the octant. That is indeed what Figs. 
1^ and 1^ show. Note also that J is proportional to the 
imaginary part of hi 2 I 123 ^ 13 ) where h = MM^ [TD]. For 
a negligible smallest mass h has a dominating 23-block 
in the normal mass ordering, whereas for the inverted 
ordering it has a democratic structure. Taking into 
account that predictions in the inverted ordering are 
in general less robust motivates to assume that the 
probability of sin 6 ^ — sin 6 is initially much larger 
than for the case of normal ordering. Indeed, see Fig. 
one starts with almost 30 % for small masses. Again, for 
quasi-degenerate masses the sign is essentially random. 

Interestingly, it is possible to change the mass ordering 
when perturbations are applied. This requires obviously 
quasi-degenerate neutrino masses, and Fig. [ 7 ] shows that 
for values around 0.1 eV the ordering can change, quickly 
reaching a probability of almost 50%. For an inverted 
ordering the probability is larger and starts for smaller 
neutrino masses. This can be traced to the larger fine- 
tuning of neutrino masses in the inverted ordering: for a 
smallest neutrino mass of 0.2 eV, we have (mi, m 2 , m 3 ) = 
(0.2,0.200187,0.206004) eV in the normal ordering, but 
(m 3 ,m 2 ,mi) = (0.2,0.205822,0.206004) eV in the in¬ 
verted one (choosing the best-values of the mass-squared 
differences). Therefore, the heaviest and next-to-heaviest 
masses are much closer together in the inverted ordering. 
After adding a perturbation, switching from inverted to 
normal is thus more likely than the other way around. 

III. PERTURBATIONS ON AN A4 MODEL 

In this section we will see how realistic our findings 
from the general case treated so far are. We apply cor¬ 
rections to a specific flavor symmetry model. 

A. The model 

We consider a model based on the discrete group A 4 , 
as developed in [IIHI3]. In the unperturbed limit, the 
charged leptons are diagonal and the neutrino mass ma¬ 
trix is 



The mass-dimension parameters b and c are related to 
vacuum expectation values (VEVs) of A 4 singlets and 
respectively. An A 4 triplet field (p' acquires VEVs in 
the ( 1 , 1 , 1 ) direction and governs the parameter d: 

{n=Ub, {0=Uc, {p')=v'{1,1,1), (4) 

with b = UbXb,c = UcXc,d = v'xd, and X},^c,d are dimen¬ 
sionless parameters. 
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Taking b, c and d as free parameters, the zeroth order 
mass matrix Eq. ^ can fit current neutrino data very 
well. To obtain the required parameters and to facilitate 
the perturbation, we minimize the following y^-function: 


X^{b,c, d) = 


023-0°23V , ' 


o’la 

m - TOsm 

Cm,™ 0 


0'23 


<^12 
2 


( Am^ - Awg ^ 




V 


^ Am^ 


(5) 


where the last term is added to fix the smallest mass, 
msni) which is implemented by taking —t 0. In 

practice, we take = Wsm/lOOO. The parameters b, 

c and d are, in general, complex numbers. We can remove 
an overall phase so only five degrees of freedom remain. 
The above mass matrix is (partly) form-invariant, the 
eigenvector to the eigenvalue 6 -I- c is always (1,1,1)^, 
hence |Cfe2p = independent of the magnitude of the 
mass matrix entries. Corrections will destroy this feature. 

As usual in models of this kind, the parameters b, c 
and d have to be somewhat tuned to get the experimen¬ 
tal values of Sm^ and Am^, which makes it technically 
difficult to find the minimum of the y^-function. The 
X^-fit gives the following conclusions: 


1. The minimal value is non-zero, Xmin ~ which 
implies reasonable agreement with current data. 
The global minima are not unique but discrete, we 
find that there are four degenerate minima with 

y2 , =3 7- 

Amm ’ 


in particular we have 0 sin^ 023 = 0.437. Due to 
the constraint |C/q, 2P = 5, one has '/ 2 \Ue 3 \ cosS ~ 
1/tan2023, leading to 5 = ±55.3°; 

3. The degeneracy between the four minima corre¬ 
sponds to (5,—0:1,2) t-)- (<5,01,2) <->■ —(<5,01,2) t-)- 
(—<5,01,2). Here <5 is the Dirac phase and 01,2 
are Majorana phases. For (<5,01,2) ^ —(<5,01,2) 
the degeneracy is obvious since it means conjugat¬ 
ing Mq. Henceforth, we name the four solutions 
(+,+), (±,-), (-,+) and (-,-) if the signs of 
(sin<5, sinoi) are (±, ±), (±,-), (-,±) and (-,-), 
respectively. The four discrete minima imply that 
for a fixed smallest mass, the model predicts defi¬ 
nite CP phases, both Dirac and Majorana. Two of 
the four solutions have positive sin^ and two neg¬ 
ative; 

4. We do not need to perturb these four zeroth or¬ 
der mass matrices separately, as the (±, ±) case is 
identical to the (—,—) case, and the (±, —) case 
identical to (—, ±); 

5. For msrn ~ 0.03 eV (normal ordering) or 0.02 eV 
(inverted ordering), there is no solution as Xmin 
creases rapidly and soon gets out of the 5cr range. 
For values below, Xmin always 3.7. The reason for 
this is the neutrino mass sum-rule TO3 — mi = m2 
(here fhi are complex masses, thus including the 
Majorana phases), which implies [T3] the relations 
mi ^ -s/Am^/y^ and m3 ^ V Am^/2, respectively. 


B. Perturbation 


2. The reason why we cannot have arbitrarily small 
X^ is because of |C/e2p = |, which forces sin^ 0i2 to 
values larger than i (sin^ 0i2 = 0.341, to be pre¬ 
cise), while the Icr-range from global fits is below 
I. The other oscillation parameters can be repro¬ 
duced to their best-fit values at the x^-minimum. 


For simplicity we study the following simple VEV mis¬ 
alignment: 

(C") = u&(l ± ei), (^') = Uc{l + € 2 ), (6) 

{ip') = v'{l ± £3,1 ± £4,1 ± £5). (7) 

As a result, the mass matrix is 


/ |d(l±£3) 6 (1 ± £ 1 ) — |<i (1 ± £5) c (1 ± £2) — (1 ± £ 4 ) \ 

A7 = I 6 (1 ± £ 1 ) — |<i (1 ± £ 5 ) c (1 ± £ 2 ) ± 3 <i (1 ± £ 4 ) - 5 *^ (1 “f £ 3 ) 1 • ( 8 ) 

yc(l±£2)— 3d(l±£4) —|(i(l±£3) 6(l±£i)±|d(l±£5)y 


Recall that the original b, c, d are fixed by our initial x^ 
minimization, so the various ei parameters are indeed 
required. Similar to our study of general perturbations, 
we randomly generate the e and study the robustness of 
the mass matrix, i.e. the stability of the octant of 623 , 
the sign of sin <5 and the mass ordering. Figs. and 
[To] show the result. To illustrate our findings, we use 
^ \e,\ < 0.04 for Fig. [S] and Fig. but X) IdI < 0.2 for 


Fig. in addition we do not give results for all zeroth 
order solutions corresponding to the signs of the phases. 

As one would expect from the general analysis in Sec¬ 
tion |IT| all percentages in the figures increase for increas¬ 
ing smallest mass. The normal mass ordering is some¬ 
what more tuned than the inverted one, as the neutrino 
mass sum-rule 7713 — rhi = m2, which holds also after 
perturbation to good precision, requires more tuned Ma- 
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smallest mass/eV 

Figure 8. Change of octant: same as Fig. for the A 4 model 
with VEV misalignment. Blue is the normal mass ordering, 
red inverted. The perturbation is made on the discrete solu¬ 
tions (-I-,-f) and ( — 



smallest mass/eV 

Figure 9. Change of the sign of sin 5: same as Fig. for 
the A4 model with VEV misalignment. Blue Is the normal 
mass ordering, red inverted. The perturbation is made on the 
discrete solutions (-1-,—) and 


jorana phases in the normal ordering M- Therefore, the 
difference between normal and inverted ordering is not 
as large as in the general case, but the overall structure 
of the plots in Figs. EE and is the same as in the 
general case. 



Figure 10. Change of mass ordering: same as Fig. for the 
A4 model with VEV misalignment. Blue is the normal mass 
ordering, red the inverted. The perturbation is made on the 
discrete solutions (-1-,—) and 

IV. CONCLUSION 

We have studied the robustness of neutrino mass ma¬ 
trix predictions in the general case and within a specific 
flavor symmetry model. We illustrate the need to include 
corrections to a mass matrix by showing that the octant 
of 623 , the sign of sin:5, or even the mass ordering can 
change when perturbations are added. Most of the re¬ 
sults are intuitively clear: 0 i 2 and <5 are more unstable 
than 013 and 623 , thus putting doubt on the discriminat¬ 
ing power of the solar neutrino mixing angle and the CP 
phase when corrections are ignored. Predictions from 
an inverted mass ordering are more unstable than the 
normal one. The larger neutrino masses are, the more 
unstable are the predictions. Going beyond 0.1 eV can 
even change the mass ordering from normal to inverted, 
quickly reaching a probability of 50 %. 

We have made here conservative assumptions about 
the perturbation parameters, namely multiplicative cor¬ 
rections. Additive corrections are expected to give qual¬ 
itatively similar perturbations, but at least as sizable as 
the multiplicative ones under study here, as they influ¬ 
ence small entries of the mass matrix more significantly. 
We have also not considered often and extensively studied 
charged lepton corrections, which are model-dependent 
and only influence the mixing matrix, independent of 
neutrino mass and mass ordering. At the current stage, 
we feel that our results already illustrate potential issues 
with the discriminative power of mass matrix predictions 
when perturbations are ignored, but at least illustrate 
quantitatively the potential impact they can have. 
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